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Abstract 

Using a code based on the Lattice Boltzmann Equation, we have 
performed numerical simulations of a turbulent shear flow. We in- 
vestigate the scaling behaviour of the structure functions in presence 
of anisotropic homogeneous turbulence, and we show that although 
Extended Self Similarity does not hold when strong shear effects are 
present, a more generalized scaling law can still be defined. 
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In the last few years there has been a growing attention on the scahng 
properties of fully developed turbulence and, in particular, on the character- 
ization of the Probability Distribution Function of the velocity increments 
6rV = Vx{x+r)—Vx{x), i.e. the velocity difference in the x— direction between 
two points at distance r. 

To this aim, usually one considers the scaling properties of the structure 
functions defined as: 

Fnir) = {\5rvn. (1) 

According to the Kolmogorov theory|l[] a scaling law for (|l]) is expected to 
hold in the so called inertial range, r/ <^ r ^ L,( L being the integral scale 
of the flow and 77 the Kolmogorov scale): 

F„(r) = A„(er)3 (2) 

where An are dimensionless constants and e is the mean rate of energy dissi- 
pation. 

There has been many experimental and numerical results suggesting that, 
because of the intermittency of the velocity field, the relation (0) is violated, 
giving an anomalous scaling law with scaling exponents Cn 7^ f • 
By taking into account the fluctuations of the energy dissipation field, the 
equation has been modified by Kolmogorov who introduced the Refined 
Similarity Hypothesis (RSH): 

S n 

F,(r) = A„(e.^)r3 (3) 
where e,. is the local rate of energy transfer. 




At present, most of the efforts, both theoretical and experimental, are de- 
voted to the determination of the anomalous scaling exponents and to the 
investigation of the role played by the RSH. 

The aim of this work is to investigate the scaling properties of the structure 
functions in the case of a homogeneous shear flow, as a simple example of 
anisotropic homogeneous turbulence. We are mainly interested to study the 
scaling laws of the structure functions and to establish if the Extended Self 
Similarity, recently introduced in literature|Q, still holds for shear flows. 
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i.e. in presence of a non isotropic turbulent flow. 

In this Letter we first remind some concepts about Extended Self Similarity 
(E.S.S.) and its relevance in order to estimate the Cn- 

Next we briefly describe the shear flows and some of their properties. Finally 
we discuss the numerical simulation and show that E.S.S docs not hold for 
shear flows, while a generalized scaling law, involving both E.S.S. and R.S.H. 
is valid. 

In principle, we can determine the scaling exponents by means of ex- 
perimental and numerical measures, but in the latter case some technical 
problems arise. 

The highest Reynolds numbers that can be achieved by laboratory experi- 
ments are about 10^"^'', while the numerical simulations performed with the 
most powerful computers now avalaible can reach i?e ~ 10^. As the computa- 
tional effort grows like Re^, it could seem very hard to obtain good estimates, 
at least comparable to the experimental results, of the scaling exponents by 
the numerical simulations. 

The concept of Extended Self Similarity (E.S.S.) can help us to fill up this 



The idea is to investigate the scaling behaviour of one structure function 
against the other, namely 



In particular it is expected that, at least in the incrtial range, f3{n, 3) = (n- 
Actually, there is strong evidence that E.S.S. is a powerful tool to investigate 
the scaling laws and that it has many advantages respect to the usual scaling 
against r, namely: 

- it holds down to the dissipative range r ~ 4 5r], 

- it holds also for low Reynolds numbers. 

Last but not least, the two previous properties allow a very accurate deter- 
mination of the scaling exponents. Indeed, the Cn can be estimated with an 
error of just a few percent. 

The above statements can be summarized as follows. We can always write 
the structure functions in the following way: 



gap. 




(4) 
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with Uq = F^lr), L = U^/e being the integral scale, and Cp dimensionless 

constants selected in such a way that fpir/rj) = 1 for r 3> r/. 

E.S.S. implies that, for all the orders p, the function fp{r /rj) = fir/rf) is the 

same. 

We want to understand which are the effects of the lack of isotropy on the 
anomalous scaling law defined in (^). To this effect, we consider a simple 
shear flow. 

Let us consider the usual Navier-Stokes equations describing a viscous, in- 
compressible fluid of density p, and velocity field v{x,t): 

dtv + {v-V)v = -^Vp + uAv + f (6) 
V ■^/= 

— » 

Let us indicate the stationary solution of the above equations as U, and 
define the turbulent velocities w as: 

v = U + w. (7) 

In order to simplify the following discussion we choose the x— direction as 
the direction of the main flow: = U, Uy = 0, = 0. 
We have a homogeneous shear flow when the main motion has a constant 
velocity in a given direction and a constant lateral velocity gradient through- 
out the whole field, e.g. Ux = U (z) and ^ = S, so there is an evident lack 
of isotropy in the system. 

Moreover we have a non zero turbulence shear stresses tensor, the compo- 
nent (wxWz) is different from zero and it makes a positive contribution only 
to dt{wl), resulting in non isotropy. 

A generalization of the "|" Kolmogorov equation for anisotropic homoge- 
neous shear flow |^ suggest that the typical scale fixed by the shear intensity 
IS Ts ~ {e/S^)2. With zero shear this scale is infinite, otherwise it has a finite 
value: below this scale the shear effects are expected to become negligible. 
The particular question we want to address is: what does it happen to the 
scaling laws (§) when falls into the inertial range? 

In order to answer this question we perform a direct numerical simulation of 
a turbulent shear flow, using a code based on the Lattice Boltzmann Equa- 
tion, for computational details see, for instance, p, |[. 
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Wc simulate a 3D fluid occupying a volume oi V = sites with L = 160, 
viscosity u = 0.014, and obeying to the usual N-S equations plus a forcing 
term / = {fx{z),0,0) chosen such that the stationary solution of the N-S 
equations is: 

= Asin{k^z) Uy^O t/^ = 0. (8) 

kz = ^ being the wave vector corresponding to the integral scales, and 
A = 0.3. 

In this way the shear has a spatial dependence S{z) ~ cos{kzz). We can 
access both zones where the shear is maximum and locally homogeneous, 
and zones where the shear is minimum. 

Wc evaluated Vrms as the mean value of The simulations have been 

done at Rex = ^^"J^" ~ 40, with A ~ 15 lattice spacings, and the Kolmogorov 
scale is about 1 lattice spacing wide. 

The simulation has advanced 100000 iterations corresponding to about 25 
macroscale eddy turnover times Tq ~ L/vrms' 40 velocity configurations have 
been saved every 2500 time steps, in order to ensure the statistical indepen- 
dence of the different configurations. 

We have evaluated the structure functions -F„(r) up to the tenth order. The 
mean values of \6rv\" have been evaluated through time and spatial average 
at fixed z-level: 

{0{r,t)) ^^£dtj^l dxdy 0{r,t). 

In Fig.l we have a log-log plot of the longitudinal (x-direction) structure 
function Fe(r) against F^^r), obtained from the velocity fields corresponding 
to the minimum shear level. The dashed curve is the best fit done in the 
range between the 20-th and 30-th grid point, and corresponds to a slope of 
1.79 in good agreement with other measured values of (q. Every point in the 
plot corresponds to a grid point and the lattice spacing is ~ Ir] wide. As we 
can see the E.S.S. holds as usual until 4 577. 

Fig. 2 shows the same plot but at the maximum shear level. It is quite evident 
that E.S.S. docs not hold. In any case, the slope corresponding to the best 
fit can be estimated at about 1.43, quite different from the previous value. 
Similar results have been obtained for all the others structure functions. 

In Table 1 we show the scaling exponents obtained for the even order 
structure functions. 
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C2 


C4 


Ce 


Cs 


ClO 




min sh 


0.70 


1.28 


1.79 


2.25 


2.68 




max sh 


0.76 


1.18 


1.43 


1.56 


1.61 




SL mod. 


0.696 


1.279 


1.778 


2.211 


2.593 



Table 1: Scaling exponents evaluated at the minimum shear (first line), at 
the maximum shear (second line), and from the She-Leveque ||1 1|| model. 



We can suggest the following explanation for the different scaling behaviour 
in presence of shear. In our simulations the scale is about 4 lattice spacings 
at the maximum shear level, so the entire range over which the E.S.S. holds 
(see Fig.l) is subjected to the shear effects. 

Our result clearly shows that the shear completely destroys the E.S.S. 
We now turn our attention to RSH. Following ^ we can consider the gen- 
eralization of RSH by introducing an effective scale S{r) = {SrV^) / (er) = 
rf{r/ri). Then ESS combined with RSH suggests: 

If the equation (^ is true, as it has already been verified for experimental 
data sets referring to homogeneous and isotropic turbulence 0, we expect 
that the local rate of energy transfer and the structure functions satisfy the 
following scaling law: 

{6^) ~ {e'^){6rvY (10) 
over a range wider than the inertial one. 

Using the data from our simulation, we obtained the results shown in Fig.s 
3-4. 

As we can see the scaling of {e"^) {ISrvf)"^ against {SrV^) is well verified in both 
the zones of maximum and minimum shear with a slope very close to one. 
This result is extremely interesting and suggests that the scaling law ( p!0D is 
universal, regardless the isotropy conditions of the turbulent flow. 
Let us summarize the results that have been obtained and suggest a possible 
interpretation for them and what should be their future developments. 
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First of all it has been shown that E.S.S. does not hold for anisotropic turbu- 
lent flows, according to similar results obtained from experimental data sets 
of turbulent boundary layers |T^, where strong shear effects are expected to 
appear. 

It means that moments of different order show a different dependence from 
the cutoff scale. This means that the shear affects the function fp{r/ri), de- 
fined in (1^), which is no longer the same for all the orders p. 
Nevertheless, the scaling law ([T0|) is valid even in presence of shear and at 
the smallest scales investigated, suggesting that the scaling law of a generic 
structure function is related to those of the third one and of the energy dis- 
sipation in a universal way, for all analyzed scales, a remarkably non-trivial 
result. 

We think that the investigation of the self-scaling properties of the energy 
dissipation would deserve more attention, in order to understand how the 
structure functions of the velocity increments depend on the resolution scale 
and to explain the ESS violation in shear flows. 

A deeper analysis of these arguments, together with other numerical and ex- 



perimental results, will be the subject for further investigation |13, 14 
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Figure Captions. 



Figure 1. Log- log plot of F^i^r) against F^lr) at the minimum shear. The 
dashed line is the best fit with slope 1.79. Every point in the plot 
corresponds to a grid point and the lattice spacing is ~ Ir] wide. 

Figure 2. The same as in Fig.l at the maximum shear. The dashed line is 
the best fit with slope 1.43. 

Figure 3. Plot of {e^) {\Srv\^)'^ against (SrV^) at the minimum shear. The 

points refer to the scales at 2,4,5,8,10,16,20,32,40 grid points and 
the dashed line is the best fit done over these points, corresponding to 
the slope 0.99. 

Figure 4. The same as in Fig.3 at the maximum shear. The dashed line is 
the best fit with slope 0.99. 
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